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We have studied the vibrational properties of CO adsorbed on platinum and platinum- 
ruthenium surfaces using density-functional perturbation theory within the Perdew-Burke- 
Ernzerhof generalized-gradient approximation. The calculated C-0 stretching frequencies are 
found to be in excellent agreement with spectroscopic measurements. The frequency shifts that 
take place when the surface is covered with ruthenium monolayers are also correctly predicted. 
This agreement for both shifts and absolute vibrational frequencies is made more remarkable by 
the frequent failure of local and semilocal exchange-correlation functionals in predicting the sta- 
bility of the different adsorption sites for CO on transition-metal surfaces. We have investigated 
the chemical origin of the C-0 frequency shifts introducing an orbital-resolved analysis of the 
force and frequency density of states, and assessed the effect of donation and backdonation on 
the CO vibrational frequency using a GGA + molecular U approach. These findings rationalize 
and establish the accuracy of density-functional calculations in predicting absolute vibrational 
frequencies, notwithstanding the failure in determining relative adsorption energies, in the strong 
chemisorption regime. 

I. INTRODUCTION 

Fuel cells are energy conversion systems of potentially high environmental benefit (l| that provide 
electricity and heat by catalytic conversion of a fuel, such as hydrogen or methanol. Despite their 
advantages, several technological obstacles have hindered the deployment of fuel-cell systems. For low- 
temperature fuel cells that use platinum as electrode material, one major limitation is CO poisoning, 
whereby CO occupies active sites on the platinum catalyst and prevents fuel oxidation Typically, 
in polymer electrolyte membrane fuel cells (PEMFCs), CO concentrations must be brought below 10- 
50 ppm to maintain an acceptable catalytic performance. For comparison, CO concentrations are 
generally on the order of thousands of ppm in reformed hydrogen fuels j^, 0|. CO poisoning is even 
more problematic for direct methanol fuel cells (DMFCs) since CO is always present in critical amounts 
as an intermediate in methanol oxidation fs^. 

Ruthenium islands on platinum catalysts have been shown to considerably attenuate CO poisoning 
0. 0, S], although the microscopic details of this phenomenon are not completely understood. Two 
main mechanisms have been proposed to explain this improved tolerance to CO. Within the bifunctional 
mechanism model, adsorbed OH species generated by water dissociation at the platinum/ruthenium edge 
promote the oxidation of CO {the promotion effect) 111! ]. According to an alternative 



view, ruthenium modifies the electronic structure of neighboring platinum atoms, reducing their affinity 
for CO {the ligand/intrinsic effect) 0, To investigate further these mechanisms of central interest 
to fuel-cell technology, it is necessary to elucidate the nature of the chemical interaction between CO 
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and bimetallic surfaces. 

In most cases, density-functional theory provides a reliable description of molecular adsorption and 
dissociation on transition metals 
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171]. However, CO adsorption on transition- 



metal surfaces is unexpectedly problematic. Indeed, at low CO coverage, local and generalized- 
gradient density-functional calculations predict CO adsorption on Pt(lll) to take place at the fee 
site, contradicting low-temperature experiments, which unambiguously indicate atop adsorption. This 
well-known qualitative discrepancy (the "CO/Pt(lll) puzzle") [18|] precludes an accurate descrip- 
tion of important phenomena, such as the surface diffusion of CO adsorbates and the thermal pop- 
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2CL and a wide body of literature exists on the subject 



In this work, we highlight and rationalize the accuracy of density-functional calculations in predicting 
the stretching frequencies of CO adsorbed on platinum and platinum-ruthenium surfaces, notwithstand- 
ing the failure in predicting the most stable adsorption site. We first present density-functional theory 
and density-functional perturbation theory results for the energetic, structural and vibrational properties 
of adsorbed CO. Second, we introduce a novel orbital-resolved force analysis to clarify the electronic 
origins of the C-0 frequency shifts as a function of the adsorption site. Last, we rationalize the accuracy 
of the stretching-frequency predictions by analyzing the influence of donation and backdonation using 
a GGA -I- molecular U model recently introduced by Kresse, Gil, and Sautet {211 1 . 



II. THEORETICAL BASIS 

The (111) transition-metal surface is modeled using a periodically repeated slab composed of four 
layers, each layer containing four atoms per supercell. A \/3 x 2 adsorption structure corresponding to 
a coverage of 1/4 of the monolayer (ML) is adopted for the CO overlayer. Atomic cores are represented 
by ultrasoft pseudopotentials [33]. The exchange-correlation energy is calculated within the Perdew- 
Burke-Ernzerhof generalized-gradient approximation (PBE-GGA) [34]. The size of the vacuum region 
separating the periodic slabs is ~13 A. We use a shifted 4x4x1 mesh with cold-smearing occupations 



351 ] (smearing temperature of 0.4 eV) to sample the Brillouin zone. Energy cutoffs of 24 and 192 Ry 
are applied to the plane-wave expansions of the wavefunctions and charge density, respectively. As 
discussed in Ref. [18], the system is not spin-polarized. Using the above slab thickness and calculation 
parameters, we verify that the adsorption energies are converged within less than 10 meV and the 
atomic forces within a few meV/A. 

The bond length and stretching frequency of CO in the gas phase are calculated to be 1.140 A and 
2140 cm^-*^ (experimental values are 1.128 A and 2170 cm"-*^). The PBE-GGA lattice parameter and 
bulk modulus of platinum are 3.993 A and 2.36 Mbar, in good agreement with experimental values of 
3.923 A and 2.30 Mbar ^]. All our calculations use fully relaxed configurations. 
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TABLE I: Adsorption energies, struetural properties, and vibrational frequeneies ealeulated using density-funetional 
theory and density-funetional perturbation theory for CO adsorbed on elean Pt(lll) surfaees. 
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TABLE M: Adsorption energies, struetural properties, and vibrational frequeneies ealeulated using density- 
functional theory and density-funetional perturbation theory for CO adsorbed on ruthenium-eovered Pt(lll) 
surfaees. 
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III. RESULTS 

We report the results of our density-functional calculations in Tables[i]and Qj] For platinum surfaces, 
the calculated atop binding energy i?ads(atop) = 1.61 eV is consistent with that reported in Ref. 
(4]]] (1.55 eV in the same adsorption structure). As a matter of comparison, the experimental heat of 
adsorption at 1/4 ML CO is 1.30 eV. The relative adsorption energy i?ads(atop) — ^ads(fcc) is calculated 
to be 0.13 eV, in accordance with the gradient-corrected relative adsorption energies (ranging from 
—0.10 to —0.25 eV) reported in Ref. [18]. As expected, our density-functional calculations favor CO 
adsorption at the threefold fee and hep adsorption sites for platinum and platinum-ruthenium surfaces, 
confirming the aforementioned disagreement with experiments. (Note that bridge adsorption of CO on 
platinum-ruthenium surfaces is predicted to be energetically unstable.) Despite this noteworthy failure, 
the bond length d{C-0) is calculated to be 1.153 A at the atop site and 1.177 A at the bridge site on 
platinum, in good agreement with experimental bond lengths (1.15ib0.05 A at both the atop and fee 
sites). Similarly, the distance d(Pt-C) from the carbon to its nearest platinum neighbor, calculated to 
be 1.864 A at the atop site and 2.029 A at the fee site, is always within experimental error (experimental 
bond lengths are 1.85±0.1 A and 2.08±0.07 A at the atop and fee sites, respectively). Note that both 
bond lengths increase with coordination. 

The full phonon spectra for CO adsorbed at the atop, bridge, hep, and fee sites on platinum and 
platinum-ruthenium surfaces are calculated using density-functional perturbation theory (DFPT) [4^. 
Within this approach, the full dynamical matrix of the system is computed exactly by solving the self- 
consistent linear-response problem describing the electron response to atomic perturbations of arbitrary 
wavelength. The DFPT spectra reported in Tables[i]and[]J]exhibit some common and expected features. 
The highest vibrational frequency in the range [1700 cm^-^, 2100 cm~^] corresponds to the localized 
C-0 stretching mode. The second highest frequency z^(M-C) in the range [300 cm^-^, 600 cm^-^] is 
related to the stretching of the metal-carbon bond. This mode is followed by two CO bending modes 
with frequencies lying 20-200 cm^-^ below v{M-C). All the other modes involving displacements of the 
heavy metal atoms are found in the frequency range [0 cm^-^, 300 cm^^]. 

We now focus on the dependence of the C-0 stretching frequency as a function of the adsorption 
site. Upon atop adsorption on platinum, the predicted v{C-0) is reduced from 2140 cm^-^ to 2050 
cm^-^, corresponding to a red shift Az^(C-O) of —90 cm^^. For comparison, the experimental stretching 



frequency, as obtained by means of sum-frequency generation (SFG) spectroscopy [39||, decreases from 
2170 cm"-^ to 2070 cm^-^, corresponding to Av{C-0) = —100 cm~^. The frequency shifts are even 
more marked at high-coordination sites: v{C-0) is predicted to be 1845 cm^-^, 1752 cm^-^, and 1743 
cm^-^ at the bridge, hep, and fee sites, corresponding to red shifts of up to —397 cm^^. These DFPT 
stretching frequencies show remarkable agreement with their SFG counterparts: z^(C-O) = 1830 cm^-*^ 
at the twofold bridge site, i^(C-O) = 1760 cm-^ at the threefold hep and fee sites, corresponding to a 
maximum red shift of —410 em^-^. Accurate DFPT frequencies are also obtained for CO adsorbed on 
platinum-ruthenium bimetallic surfaces. Indeed, the calculated stretching frequencies 1979 cm^-^ (1 Ru 
ML) and 1969 cm^-^ (2 Ru ML) at the atop site compare very closely to the SFG result of 1970 em^^. 

In conclusion, all calculated CO stretching frequencies deviate by less than 2% from the measured 
ones, irrespective of the adsorption site and nature of the metal surface. The correct prediction of 
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TABLE III: Moleeular oeeupations and bond order for CO in the gas phase and for CO adsorbed on platinum. 

the frequency red shifts allows the direct recognition of CO adsorption sites and confirms that CO 
preferentially occupies atop sites on platinum-ruthenium bimetallic surfaces. This very close agreement 
with experiment is made more remarkable by the lack of accuracy of the PBE-GGA adsorption energies. 
In the remainder of this work, we show how this accuracy can be rationalized in terms of the hybridization 
of the CO molecular orbitals with the metal bands. 



IV. DISCUSSION 
A. Electronic Origins of the Frequency Shifts 

The hybridization of the metal d bands with the Itt* lowest unoccupied molecular orbitals (LUMOs) 
and the ba highest occupied molecular orbital (HOMO) plays a predominant role in the adsorption 
energy of CO on transition metals. According to the BIyholder model [4^], these electronic interactions 
result in electron donation (i.e., partial depletion of the 5a orbital) and electron backdonation (i.e., 
partial filling of the 27r* orbital). It has been shown that the trends of the adsorption energies of CO 
on transition-metal surfaces can be correlated to the amounts of donation and backdonation (Hammer- 
Morikawa-N0rskov model) [12]. Among the suggested solutions to the "CO/Pt(lll) puzzle" — e.g., 
nonequivalent GGA description of different bond orders incorrect singlet-triplet CO excitation 



energies 



23l ]. effect of metal semicore polarization [24|] — Kresse, Gil, and Sautet have proposed that 
the inaccuracy of density-functional calculations in determining the most stable adsorption site is due 
to an overestimation of the interaction between the lir* orbitals and the metal bands, resulting from 



an underestimation of the HOMO-LUMO gap [2l||. As discussed in the next section, this interpretation 
recovers the essential features of CO adsorption on transition metals: it identifies the tendency of local 
and generalized-gradient DFT to delocalize and overhybridize electronic states. Nevertheless, as shown 
below, the site-dependence of the C-0 bond length and vibrational frequency is not affected by the 
LUMO and HOMO hybridizations. In other words, the hybridizations that subtly determine the relative 
CO adsorption energies do not influence the structural and vibrational predictions. 

To establish this fact, we first introduce a spectral force analysis. The main objective of this analysis 
is to separate and assess the contribution from each CO molecular orbital to the force F acting on a 
given atom. The central quantity we introduce is 0/ i^-(e), the force density of states (FDOS) of the 
orbital x< which is defined as the x-fesolved density of states weighted by the wavefunction contribution 



6 



^a) Atop DOS 




20 

9 (eV-^; 

(b) Fee DOS 

5 



12 12 

9x (eV-i) (eV-^) 




12 
9x (eV-i) 



1 2 
9x (eV-i) 



FIG. 1: Total density of states, density of states projected on the Pt atomic orbitals, and density of states 
projected on the CO molecular orbitals for atop and fee adsorption of CO on Pt(lll). 
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FIG. 2: Total force density of states, force density of states projected on the Pt atomic orbltals, and force density 
of states projected on the CO molecular orbltals for atop and fee adsorption of CO on Pt(lll). 
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FIG. 3: Correlation between bond length and intramolecular frequency for CO on clean and ruthenium-covered 
Pt(lll) surfaces. 



to the force F/ acting on atom /. To be more explicit, the FDOS of a given CO molecular orbital x 
can be expressed as: 

0,,^(e)=^F,,,Kx|V«)l'5(e-e^), (1) 

i 

where denotes the electronic wavefunction, is the electronic energy, and 

F/_j = —fif\ipi\^dv/dYij is the wavefunction contribution to the force Fj. (The calculation 
of the overlap (xlV'i) in the ultrasoft formalism is detailed in AppendixO) It should be noted that, by 
summing the integrated FDOS 4>/^;^(e) over a complete set of orbitals satisfying orthonormality, one 
obtains the total electronic force acting on atom /. As a consequence, the FDOS can be quantitatively 
connected to relevant observables. Additionally, by projecting the force density of states along the 
normalized atomic displacements A^~^ corresponding to the C-0 stretching mode, we obtain the force 
density of states along the stretching mode = ^/ A^~^ ■ 0/^^(e), to be heuristically identified 

as the orbital contribution to the intramolecular force. 

The orbital-projected density of states (DOS) ^^^(e) = X^i KxlV'i)P<^(e — e^) is commonly used to 
provide an insightful picture of the electronic hybridizations that take place when CO is adsorbed on 
platinum. Similarly, "l^^~^(e) describes the influence of electronic hybridizations on the force along the 
C-0 stretching mode. The projected densities of states Qxi^) and projected force densities of states 
<t>^~^(e) for different adsorption sites are plotted in Figures [Hand [21 A detailed analysis of the orbital- 
resolved densities of states is given in Ref. j2l|. For the purpose of our study, we emphasize the 
following features. When CO adsorbs on Pt(lll), the 4a* and 5a orbitals hybridize with the metal d^2 
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band, generating 4a* and 5a states with mainly adsorbate character (adopting the terminology of Ref. 
(44I . the tilde symbol denotes hybrid states). These 4a* and 5a states are found in the energy ranges 
[—12 eV, —9 eV] and [—9 eV, —5 eV] relative to the Fermi level. Above —5 eV, the 5a orbital and the 
d^2 band generate a d^r band with predominant metal character. This interaction results in a partial 
depletion of the 5a HOMO (electron donation). In addition, the interaction between the Ivr and 27r* 
orbitals and the dxz and dyz bands produces Itt states in the range [—9 eV, —5 eV] and a broad d^ 
band above —5 eV, causing partial occupation of the Itt* LUMOs (electron backdonation) (43!]. The 
changes in molecular orbital occupations due to CO adsorption are reported in Table [TTTl 

In order to understand how the generation of these hybrid states affects the intramolecular force, we 
turn to the FDOS (Figure [2]). The graphs are plotted according to the convention that bonding states 
(i.e., opposed to the stretching of the C-0 bond) correspond to negative values of "t"^"^. First, we 
note that the bonding contribution from the 3a state does not vary with the adsorption site, confirming 
that the 3a state retains a strong molecular character. Additionally, we observe that the 5a and Itt 
states are bonding while the 4a* is antibonding, as expected intuitively. In the energy region above —5 
eV, another contribution appears. This contribution corresponds to high-energy wavefunctions located 
inside the platinum slab, as evidenced by the absence of any molecular-orbital force contribution above 
—5 eV. Nevertheless, due to their metal character, the contribution of these high-energy wavefunctions 
is mostly canceled by the positively charged platinum cores. Consequently, the local contribution from 
the hybrid states of strong molecular character prevails. 

Besides these observations, the main feature of the FDOS graphs is the predominant bonding con- 
tribution between —9 eV and —5 eV. At the atop site, the curve displays a sharp negative peak which 
corresponds mainly to the Ivr orbital-resolved contribution 0f~°(e). At the fee site, both the magnitude 
and the relative share of the peak are reduced, clearly indicating that the Itt states have more influence 
on the change in intramolecular bonding than any of the other hybrid wavefunctions. The Itt states 
maintain a predominant Ivr character at the atop site, whereas at high-coordination sites this molecular 
character is significantly reduced due to a stronger hybridization with the substrate. Therefore, the Itt 
bonding contribution to the intramolecular force decreases with site coordination. As intramolecular 
bonding decreases, the C-0 bond length increases. The predominance of the Ivr bonding contribution 
is confirmed by the density-distribution analysis initially introduced by Zupan, Burke, Ernzerhof, and 
Perdew [451], as discussed in Appendix [Bl 

For CO adsorbed on transition metal surfaces, the intramolecular bond length and the intramolecular 
stretching frequency are strongly correlated. An extensive study of Gajdos, Eichler, and Hafner [4^ 
showed a linear correlation between d{C-0) and z/(C-0) for CO adsorbed on close-packed transition 
metals: i^(C-O) shifts down in frequency as d{C-0) increases. As illustrated in Figure O a similar 
trend is observed for CO adsorbed on ruthenium-covered platinum surfaces. Therefore, the increase in 
C-0 bond length at high-coordination sites, which reflects a decrease in Ivr bonding contribution, is 
accompanied by a reduction of the C-0 stretching frequency. 



While the preceding is consistent with the interpretation given in Refs [43, |44j, |46|| , it is important to 
make one central observation: although the LUMO 2vr* filling is a reasonable measure of the amount 
of hybridization between the Ivr, 2vr* orbitals and the metal d^z, dy^ bands, filling the 2vr* orbitals does 
not directly weaken the bond, as evidenced by the very low values of ^^^f^ in the energy range [—9 



10 



eV, —5 eV]. This interpretation helps explain the fact that the CO adsorption energies do not show a 



well-defined relationship with the C-0 stretching frequency [471. l48l] . 

The main conclusion of this section is as follows. At variance with the CO adsorption energies, 
electron backdonation and electron donation have little direct bearing on the intramolecular forces. 
Their immediate effect on the molecular bond length and stretching frequency cannot account for the 
observed shifts. Instead, the changes in bond length and stretching frequency are primarily related 
to the hybridization of the Ivr molecular orbitals. This provides important indications as to why the 
structural and vibrational properties of CO adsorbed on platinum and platinum-ruthenium surfaces are 
accurately predicted. The GGA -|- molecular U study presented in the next part provides additional 
quantitative evidence in support to this conclusion. 



B. Influence of Donation and Backdonation on the Accuracy of the Frequency Predictions 

As mentioned above, the failure of density-functional calculations in predicting CO adsorption en 



ergies is traceable to an overhybridization of the CO molecular orbitals with the metal bands [2ll. l49l|. 
To assess the influence of this overhybridization on the accuracy of the calculated adsorption energies, 
bond lengths, and vibrational properties, we have performed a sensitivity analysis. This analysis con- 
sists of controlling and varying the HOMO and LUMO hybridizations, while monitoring the variations 
of the mentioned observables. To this end, we have used the GGA -|- molecular U approach introduced 
by Kresse, Gil, and Sautet jill. This approach (inspired by the LDA -|- U method [50]) consists of 
adding an orbital-dependent term to the GGA energy functional, thus imposing a penalty on orbital 
hybridization. 

We employ the following GGA -|- molecular U energy functional: 

EgGA+U = EgGA H 1^ ^ Tr{f2,r*,o-(I - f27r*,o-)} + ^ /5(7,o-(l - /Scr.cr), 

cr cr 

where /sct,^ is the occupation of the 5cr orbital of spin a and ^21^*, a is the occupation matrix associated 
with the 2tt* and 27r* orbitals of spin a. The parameters U2tt* and U^a penalize noninteger occupations 
of the 2-K* and 5cr orbitals: U2tt* reduces 27r* backdonation while U^^ reduces 5cr donation, as illustrated 
in Figure[4l The parameters U2tt* and U^cf can also be interpreted as shifting the effective single-electron 
energies. Heuristically, modifies the single-electron energy by an amount C/^(l/2 — /^). Thus, 
C/27r* increases the 2tt* energies, whereas U^a decreases the 5cr energy, causing the HOMO-LUMO gap 
to increase. 

The present functional differs slightly from that introduced by Kresse, Gil, and Sautet. The GGA 
-|- molecular U energy in Eq El whose expression is based on the matrix formulation introduced by 
Cococcioni and de Gironcoli [Sll], is invariant with respect to the choice of the x- and y-axes. In other 
words, an arbitrary rotation of the molecular orbitals does not affect the GGA -|- molecular U energy. 
Additionally, the functional allows the freedom to vary both the amount of electron backdonation and 
that of electron donation. The necessity of simultaneously varying backdonation and donation will be 
discussed later. 

Although a molecular U term is admittedly a simplified energy correction, it reproduces the essential 
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FIG. 4: Density of states projected on the CO molecular orbltals as a function of f727r* and Usc^ controlling the 
hybridization of the LUMO and HOMO orbltals for atop adsorption of CO on Pt(lll). 
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features of the energetics of CO adsorption [21']. The cluster calculations of Gil et al., based on the 
B3LYP hybrid functional [52], confirm that the inaccuracy of the GGA energies can be ascribed to 
an overhybridization of the 2tt* molecular orbitals |49]. This conclusion is supported by the recent 
periodic-slab B3LYP calculations of Neef and Doll |i25|, |26t|. Moreover, experimental studies indicate 
that the adsorption energy of CO on platinum shows a linear dependence with respect to the energy 
of the center of the metal d bands, 35 in agreement with the theoretical model developed by Hammer, 
Morikawa and N0rskov 12]- However, the coefficient of proportionality is overestimated within density- 
functional calculations, indicating that the interaction between the 2tt* orbitals and the metal d bands 
is excessive. 

We thus proceeded to calculate the energetic, structural, and vibrational properties for CO adsorbed 
on platinum. Stretching frequencies are now obtained by diagonalizing the two-by-two dynamical matrix 
associated off-equilibrium displacements of the carbon and oxygen atoms in the direction normal to the 
surface. Due to the large atomic mass of platinum, the resulting stretching frequencies deviate by less 
than 1 cm^-^ from the full DFPT phonon frequencies. 

The results of the calculations are presented in FigureO and in Appendix[Cl along with methodolog- 
ical details. As expected, adsorption energies decrease with increasing penalization on the hybridizations 
of the HOMO and LUMOs. Consequently, both donation and backdonation favor CO adsorption, in 
agreement with Ref. [21]. Moreover, we observe that electron backdonation tends to decrease the 
relative adsorption energy £'ads(atop) — -Egdstfcc), confirming that 2tt* backdonation favors CO adsorp- 
tion at high-coordination sites, as demonstrated by Anderson and Awad [53]. Additionally, the effect 
of Usa on the relative binding energy is much weaker than that of C/27r*. This result supports the 
hypothesis that the failure of density-functional calculations in predicting the most stable adsorption 
site is principally related to an overestimation of 2tt* backdonation 0, Eq, 21|. 



Considering now the structural and vibrational properties, we observe more complex U2tt*- and Usa-- 
dependencies. The effect of electron donation must clearly be taken into account when analyzing the 
sensitivity of the calculated stretching frequencies. Note that the invariance of the stretching frequency 
with respect to U2tt* for CO adsorbed at the atop site, as already observed in Ref. [3] for copper 
surfaces, can be explained by the fact that the [/27r*-axis is tangent to the contour line i^(C-O) = 
2050 cm^-^. Despite this fact, the dependence of z^(C-0)(atop) with respect to is appreciable, 
supporting the idea that 27r* backdonation alone does not control the site-dependence of the C-0 
stretching frequency. 

To conclude this section, we report the results of our sensitivity analysis (FigureE]). Large ranges for 
[/27r*and Usa are selected: eV < U2tt* < 5 ev and —10 eV < U^a < 10 eV. Note that the parameter 
U2tt* is kept positive since the GGA + molecular U calculations clearly indicate that 27r* backdonation is 
overestimated. We emphasize that these energy ranges correspond to large shifts in the single-electron 
energies (up to Ae2n* = 1 eV and lAes^l = 4 eV) and to large variations of the adsorption energies (up 
to |AE'ads| = 1-5 eV). Thus, the relative variations of the adsorption energies are comparable to their 
absolute values. Despite these sizable variations of the adsorption energies, we observe little variations 
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FIG. 5: Adsorption energy and intramolecular stretching frequency as a function of \J2-n* and U^^r controlling 
the hybridization of the LUMO and HOMO orbitals for atop and fee adsorption of CO on Pt(lll). 
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FIG. 6: Ranges of variation (indicated by 
black error bars) of the adsorption energy, 
bond length, and stretching frequency of 
CO on Pt(lll) for a very broad range of 
hybridizations (0 eV < U2,t* < 5 eV and 
-10 eV < Usa < 10 eV). 




atop ■ fee I 



of the bond lengths and stretching frequencies: 

' 1.127 A < (i(C - 0)(atop) < 1.165 A 

1.173 A < d{C- 0)(fcc) < 1.199 A 

1933 cm-i < u{C - O)(atop) < 2055 cm'^ 

1666 cm"^ < u{C - 0)(fcc) < 1878 cm"^ 

These small variations account for the remarkable accuracy of the bond lengths and stretching 
frequencies calculated within PBE-GGA. In particular, they justify the correct ordering of the C-0 
stretching frequencies despite important qualitative errors in predicting the relative CO adsorption 
energies. These results provide strong support to the conclusion of the preceding section: the variations 
of d(C-O) and z^(C-O) are not directly due to electron donation and backdonation, but rather to the 
hybridization of the Ivr orbitals. 



V. CONCLUSION 



This study has evidenced that the PBE-GGA predictions for the stretching frequencies of CO adsorbed 
on platinum are in excellent agreement with SFG experiments despite the well-known failure of local 
and generalized-gradient calculations in predicting the most stable adsorption site. Similar agreement 
is obtained for CO adsorbed on platinum-ruthenium bimetallic surfaces, allowing the direct recognition 
of CO adsorption sites. 

Our orbital-resolved force analysis has demonstrated that the variations of bond length and stretching 
frequency as a function of the CO adsorption site are principally due to the Ivr hybridization, rather 
than the 2it* and 5a hybridizations. Using the GGA -|- molecular U approach, we have performed a 
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sensitivity analysis to quantify the influence of the 2it* and 5a hybridizations on the structural and 
vibrational properties for CO on platinum. The effect of 2tt* backdonation has been shown to be small 
and comparable to that of 5cr donation, contradicting the widespread idea that backdonation controls 
the frequency shifts. 

These results explain the remarkable accuracy of the PBE-GGA frequency predictions despite errors 
in the hybridizations of the 2it* and 5a orbitals. Furthermore, they suggest a promising way to connect 
density-functional calculations with experiments in some of the most problematic cases of molecular 
adsorption on transition metals. 
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APPENDIX A: ULTRASOFT OVERLAPS 



The Vanderbilt ultrasoft formalism 33] consists of replacing the density operator n(r) with: 

^^^(r) = n(r) + ^n7(r) (Al) 



where n/(r;ri,r2) = Y^n,mQi,n,m{r)(3}n{^i)Pi,m{r2) is the charge-augmentation contribution from 
the ionic core /. Correspondingly, the overlap operator S becomes: 

S = i + Y,Si (A2) 

where 5'/(ri,r2) = J2n,ml Qi,n,m{r)drf3}^^{ri)f3i.m{r2) is the ionic contribution to the overlap oper- 
ator. The ultrasoft pseudopotential of the ionic core / is the sum of a local part vf{r) and a nonlocal 

part vfHri,r2) = En,™ ^L„./?/*,n(ri)/3/,m(r2). 

This ultrasoft formalism considerably improves the convergence of density-functional algorithms 
with respect to the energy cutoffs applied to the plane-wave expansions of the wavefunctions and 
charge density. However, to calculate the overlap {x\S\ip) between the molecular orbital x the 
wavefunction tp, it must be borne in mind that the overlap operators and corresponding to x 
and ^ are distinct since the molecular orbital is calculated without the platinum slab. The procedure 
employed here consists in including fictitious platinum cores in the calculation of the molecular orbital. 
These fictitious cores are obtained by setting the local and nonlocal part of the platinum pseudopotential 
to zero, while keeping the contribution to the overlap operator unchanged. The resulting operator S'^*'^* 
being identical to S^, the overlap coefficient can be calculated as (x-^*^*|'5x*'^*|V') = (X'^*^*l'S'v)|V') where 
^/«ct is the molecular orbital calculated in the presence of the fictitious platinum cores. The primary 
advantage of this procedure is that it only requires changing the pseudopotentials. 
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FIG. 7: Derivative of the density distribu- 
tion function A''-~'^'(7i^x(rs) and derivative 
of the density-gradient distribution func- 
tion A('^"°)g3,^(s) along the C-0 stretch- 
ing mode for each molecular orbital x- 
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APPENDIX B: DENSITY-DISTRIBUTION ANALYSIS 

As an additional confirmation to the predominance of the Ivr contribution in the C-0 intramolecular 
force, we have calculated the density distribution and density-gradient distribution [45^] projected on 
each CO molecular orbital x- We define the projected density distribution as: 

giA^s) = E^^l^^l^^)!' / m^tKrs-rs{v))dv, (Bl) 

where is the Seitz radius. We propose a similar definition for the projected density-gradient distribu- 
tion: 

53,x(5) = E/*KxlV'*)l' / \Mr)fS{s-s{r))dr, (B2) 

where s = \Vn\/2kFn is the reduced density gradient. The derivatives of the distribution functions 
along the C-0 stretching mode A^'~~^'>gn,x = d{C-0)dgn,x/9d{C-0) (n=l,3) are plotted in Figure[71 
Considering A^^'^'i^i .^(r^), we observe that increasing d(C-O) tends to decrease the electronic 
charge in spatial regions of high electronic density. This trend is particularly marked for the Ivr molecular 
orbitals, confirming their predominance in the intramolecular force. The predominant Ivr contribution 
can also be seen in the A^'~~^^ g^ y.{s) graph. It is important to note that the Ivr orbitals result in a 
significant increase in charge-density homogeneity. According to the bond-expansion criterion derived 
by Zupan, Burke, Ernzerhof, and Perdew (Eq. 9 in Ref [45]), this observation confirms that the Ivr 
orbitals are strongly bonding. 
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APPENDIX C: GGA + MOLECULAR U FORCE DENSITY OF STATES 

The calculation of the FDOS within GGA + molecular U can be performed along the same general 
lines as the method presented in Sec. IIV Al the only modification being in the expression of the 
wavefunction contribution to the force F/ In the ultrasoft formalism and in the absence of molecular 



U contribution, the wavefunction contribution to the force can be written as 54l. bSl] 



dv dS 



Adding the molecular U term, the force contribution becomes: 



9R/ 



(C2) 



where ^2-K*,i(T = [/io-('0icr|'S'|27r*)(27r^|S'|?/;jo-)]Q^ denotes the contribution from the wavefunction ipi^j to 
the 27r* occupation matrix i2-K*,a and fsa,ia = fia{'4^ia\S\ba) {ba\S\il}ia) is the contribution from V'io- to 
/scTCT- Note that the derivatives must be calculated keeping V'io- fixed. 

The main computational difficulty in determining the derivatives of the occupation coefficients is 
the evaluation of the response of the molecular orbitals x to the atomic displacements dx/dHj. These 
responses can be obtained by performing a separate calculation for an isolated CO molecule using the 
same linear-response approach as that employed in the phonon calculation [42]. The projected and total 
FDOS spectra as a function of the hybridization parameters for CO adsorbed at the atop site on Pt(lll) 
are reported in FiguresOandO The observed trends confirm that large shifts in the LUMO and HOMO 
single-electron energies do not affect the predominance of the In contribution. For positive values of 
U2n* and Us(7, the magnitude of the peak in the Ivr FDOS increases due to the stronger molecular 
character of the Ivr orbitals. Note that this increase in the Ivr bonding contribution is accompanied 
by an increase in the antibonding contribution from the metal bands. Similar compensation effects in 
the electronic forces can be observed for negative values of the hybridization parameters. In particular, 
for U2tt* equal to —5 eV corresponding to an unphysical overestimation of the 2it* hybridizationthe 
decrease in the bonding contribution from the Ivr orbitals is partially offset by the decrease in the 
antibonding contribution from the d bands. This observation provides a more specific understanding of 
the invariance of the stretching frequency as a function of the HOMO and LUMO hybridizations. 
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FIG. 8: Force density of states projected on the CO molecular orbitals, and total force density of states as 
functions of U2-K' controlling the hybridization of the LUMO orbitals for atop adsorption of CO on Pt(lll). 
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FIG. 9: Force density of states projected on the CO molecular orbitals, and total force density of states as 
functions of Usa- controlling the hybridization of the LUMO orbitals for atop adsorption of CO on Pt(lll). 
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